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Abstract. 

A comprehensive theory of interfacial fluctuation effects occurring at 2D wedge 
(corner) filling transitions in pure (thermal disorder) and impure (random bond- 
disorder) systems is presented. Scaling theory and the explicit results of transfer matrix 
and replica trick studies of interfacial Hamiltonian models reveal that, for almost 
all examples of intermolecular forces, the critical behaviour at filling is fluctuation- 
dominated, characterised by universal critical exponents and scaling functions that 
depend only on the wandering exponent C. Within this filling fluctuation (FFL) regime, 
the critical behaviour of the mid-point interfacial height, probability distribution 
function, local compressibility and wedge free energy are identical to corresponding 
quantities predicted for the strong-fluctuation (SFL) regime for critical wetting 
transitions at planar walls. In particular the wedge free energy is related to the SFL 
regime point tension which is calculated for systems with random-bond disorder using 
the replica trick. The connection with the SFL regime for all these quantities can be 
expressed precisely in terms of special wedge covariance relations which complement 
standard scaling theory and restrict the allowed values of the critical exponents for 
both FFL flUing and SFL critical wetting. The predictions for the values of the 
exponents in the SFL regime recover earlier results based on random-walk arguments. 
The covariance of the wedge free energy leads to a new, general relation for the SFL 
regime point tension which derives the conjectured Indekeu-Robledo critical exponent 
relation and also explains the origin of the logarithmic singularity for pure systems 
known from exact Ising studies due to Abraham and co-workers. Wedge covariance is 
also used to predict the numerical values of critical exponents and position dependence 
of universal one-point functions for pure systems. 
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1. Introduction 



Fluids adsorbed near wedges, cones and corners sliow filling phenomena |1|, ^ ^, |, 
% 0' B B E' mi similar to the wetting of planar wall-fiuid and fiuid-fiuid interfaces 
|13|. Above a filling temperature Tfiu a wedge (say) in contact with vapour at 
bulk coexistence is completely filled by liquid, analogous to the complete wetting of 
a planar wall-fiuid interface above the wetting temperature T„et- There are, however 
a number of notable distinctions between wetting and the different possible types of 
filling transition. Firstly, thermodynamic arguments |^, |^ dictate that filling precedes 
wetting, occurring when the contact angle satisfies 9 = a, where a is the wedge tilt angle 
(see figure [^). Thus filling may occur in the absence of any wetting transition i.e. if the 
walls are partially wet up to the bulk critical temperature. Secondly the conditions for 
observing continuous (critical) filling 0, ^ in the laboratory are much less restrictive 
than for continuous (critical) wetting . For example 3D cone or corner filling should 
be continuous provided the line-tension is negative. Finally, the critical singularities and 
fiuctuation effects occurring at critical filling refiect the divergence of different length- 
scales compared to wetting. For example, in a 3D wedge, long-wavelength fluctuations 
in the interfacial height along the wedge dominate and lead to an interfacial roughness 
that is much larger than for wetting at a planar wall and which exhibits universal 
properties]^, |^, [1^ . 



a 



hq(-^) 
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Figure 1. Schematic portrait of the equihbrium interfacial height leq{x) for a fluid 
adsorbed in a wedge close to a filling transition, denotes the mid-point interfacial 
height. 



For 2D wedges (and 3D cones) on the other hand something quite different and 
unexpected occurs. Recent studies based on effective interfacial Hamiltonians 0, |^, ^ 
and more microscopic Ising models indicate that there is a fundamental connection 
with the strong- fluctuation (SFL) regime of critical wetting [|13], |16|, [1^. More 

speciflcally, even in the presence of long-ranged forces, the divergence of the interfacial 
height at critical fllling together with the precise scaling form of the mid-point interfacial 
height probability distribution function (PDF) are identical to predictions for the critical 
wetting SFL regime. In other words the substrate geometry effectively turns off the 
influence of long-ranged solid-fluid and fluid-fluid forces so the fluid mimics fluctuation 



Wedge covariance for 2D filling and wetting. 



3 



behaviour predicted for planar systems with short-ranged forces. In 2D this has been 
demonstrated for both pure and impure systems (with random-bond disorder) and 
can be precisely expressed in terms of special wedge covariance laws which relate the 
interfacial heights and PDF's in the different geometries at bulk two-phase coexistence 
0, |[. Whilst the demonstration of these laws, particularly for impure systems, is 
rather technical, the final expressions are extremely simple and contain a great deal of 
information. Consider a fluctuation-dominated filling transition occurring in a 2D wedge 
and let lw{9,a) and Pw{l;6,a) denote the equilibrium mid-point interfacial height and 
corresponding PDF respectively. Now let 1^(0) and Pjril', 0) denote the interfacial height 
and PDF for a SFL regime wetting transition, occurring at planar wall-fluid interface, 
written in terms of the contact angle. Covariance for the interfacial heights implies that 
at bulk coexistence and as ^ — a, 

Ue,a) = ue-a) (1) 

The statement of covariance for the PDF's is even stronger: 

P^{l-9,a) = P^{l-e~a) (2) 

implying that not only the interfacial height but also the roughness (and all moments 
of the distribution) at filling and wetting are similarly related. 

In this paper we further investigate the connection between filling and wetting and 
demonstrate covariance relations for other quantities of interest. Again, because the 
derivation of these results is rather technical we quote their final form here. For the 
excess free energy of the wedge fw{0,a) at coexistence we find that for both pure and 
impure systems 

Ue,a) = T{e)-T{e-a) (3) 

where t{6) denotes the point tension near a SFL regime critical wetting transition. 
We also consider the local compressibility for filling Xw{z',0,<y) corresponding to the 
derivative of the mid-point density profile w.r.t. chemical potential /i evaluated at two- 
phase coexistence. This is related to the corresponding expression Xn{z;9) for SFL 
regime wetting by 

Xw{z;e,a) = ( - - l] xA^',^ - Oi) (4) 



a 

These covariance relations complement the standard scaling hypothesis for the 
singular contribution to the corner free-energy and may be understood heuristically 
by considering the special influence that the wedge (and also cone) geometry has on 
interfacial configurations. The relations are rather restrictive and from them we may 
deduce the allowed values of the critical exponents for both 2D filling and critical 
wetting transitions without explicit model calculations. This is somewhat analogous 
to the restrictions imposed by the principle of conformal invariance for bulk and surface 
critical phenomena |TB|. We will also use the covariance relations to derive new results 
for the point tension and position dependence of the local compressibility. Indeed for 
the point tension we will be able to derive a conjectured critical exponent relation due to 
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Indekeu and Robledo |T^, |20| and also explain why, for pure systems, the point tension 
shows a logarithmic singularity as ^ — > [pT| . 

Our article is arranged as follows: In §11 we present the necessary background theory 
for critical wetting describing, in turn, the derivation of critical exponent relations 
and fluctuation regimes from heuristic scaling theory and the scaling of the density 
profile, probability distribution function and their short- distance expansions. In §111 we 
discuss the singularities of the point tension first recalling the hyper-scaling conjecture 
of Indekeu and Robledo fl^ and exact Ising model results for pure systems due to 



Abraham, Latremoliere and Upton |2T|. We then present a lengthy calculation of the 
point tension r for critical wetting in pure and impure systems systems using continuum 
interfacial models. As far as we are aware r has not been calculated before for systems 
with random-bond disorder and the derived expression shows singular behaviour as 
^ — >■ in full agreement with the Indekeu-Robledo conjecture. In §IV we begin our 
discussion of 2D filling transitions and first discuss critical exponent relations, heuristic 
fluctuation theory and the scaling of the density profile and PDF, which parallels our 
earlier discussion of critical wetting. In §V we recall the main results of explicit effective 
Hamiltonian studies of filling which will be further developed here. These calculations 
support the scaling theory developed in §IV and also demonstrate the covariance laws 
(I3)"(3) quoted above. The consequences of these relations and the restrictions they place 
on the allowed values of the critical exponents at 2D filling and wetting (in systems 
with short-ranged forces) are discussed at length. A new and very simple result for the 
point tension is also discussed from which we can derive the Indekeu-Robledo exponent 
relation. In §VI we turn our attention to scaling behaviour occurring off bulk coexistence 
and focus on the position dependence of the local compressibility leading to the final 
covariance relation (^). Finally we revisit the nature of filling and wetting in pure 
systems and use the covariance relations to re-derive the specific results of the transfer 
matrix studies. We finish our article with a brief summary of our main results and 
discussion of future work. 



2. 2D critical wetting and SFL regime 

We begin with a brief survey of the central results from the theory of critical wetting 
paying special attention to critical exponents, fiuctuation regimes, and the scaling of 
the density profile and PDF. Further details and original references can be found in the 
excellent review articles [O, IT 



2.1. Critical exponents and exponent relations 

Consider the interface between a planar substrate (wall), located in the z = plane, 
and a bulk fiuid phase at temperature T and chemical potential /i (and corresponding 
pressure p). The two-phase coexistence line is denoted /i = UsaxiT) with corresponding 
bulk liquid and vapour phases densities, pi and p^, respectively. Throughout this article 



Wedge covariance for 2D filling and wetting. 



5 



we set ksT = 1 for convenience. We suppose the wall-vapour interface preferentially 
adsorbs the liquid phase and is completely wet above a wetting transition temperature 
Twct- Thus the wetting phase boundary is defined by the vanishing of the contact angle 

e{T) = 0; T > T„et (5) 

The wetting transition corresponds to a singularity in the excess free-energy (surface 
tension) of the planar wall- vapour interface a^^. This is defined |^ by subtracting the 
bulk contribution from the total grand potential VL 

Cfwv = (,Dj 

in the limit of infinite volume V and planar area A. On approaching a critical wetting 
transition the adsorption T ^ {pi — Pv)Itt with 1^, the equilibrium height of the liquid- 
vapour interface above the wall, diverges continuously. The critical wetting transition 
has two relevant scaling fields which we can identify (for fixed strength of the wall-fluid 
intermolecular potential) as t' = (T„et— 7')/T„et and h = (pi — Pt,)(/isat — yu)- From scaling 
theory we anticipate that off two phase coexistence a^v contains a singular term, /sing 
which vanishes at the critical wetting phase boundary and can be written 

fsin^ = t'^-''^W^{ht'-^) (7) 

where as, A are the surface specific heat and gap exponents respectively, Wt,{x) is the 
scaling function and we have restricted our attention to t' > 0. Now at coexistence we 



have, by definition /sing = (Jwv — (o"w)/ + criv) so that from Young's equation it follows 
that the contact angle vanishes as 

9 ~ t'(2-"»)/2 (8) 

and can be used as an alternative measure of the temperature-like scaling field. In 
addition to the mean interface height 1^^ we also need consider the divergence of the r.m.s 
interfacial width or roughness C,± and the transverse correlation length measuring the 
correlations in fluctuations of the interfacial height along the wall. At h = 0, the 
divergence of these length-scales is characterised by critical exponents defined by 

and we expect scaling behaviour similar to (|^) off coexistence. Thus, along the wetting 
critical isotherm (T = T„et, /i — * 0) we define 

~ h-^ (10) 

and anticipate that the critical exponent iIj = /?s/A. Critical exponent relations 

immediately follow from the scaling hypothesis. Firstly from the Gibbs adsorption 
equation |12|, ^ 

dfsing/dh = 1^ (11) 

we have 

A = 2- as + ps (12) 
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which identifies the gap exponent. Secondly, from the compressibihty sum- rule (see for 
example |12[) 

d^fsinjdh^ oc (13) 

we have 2 — as — 2A = —2u\\ so that on eliminating the gap exponent we arrive at the 
important relation 

2 - = 2z/|| - 2/?, (14) 

which is valid for all dimensions, ranges of forces and fiuctuation regimes of interest. 
The perpendicular and transverse correlation lengths are related through the wandering 



exponent defined by |15 



e± ~ (15) 

with the value of C > dependent on the dimensionality and qualitative type of disorder. 
For discussions of wetting and filling in 2D systems the most relevant values are ( = 1/2 
and ( = 2/3 for pure (thermal disorder) and impure systems (random-bond disorder) 
respectively. Recall that random-fields destroy phase coexistence in two dimensions 
so cannot be considered. Also whilst values of the wandering exponent < 1/2 have 
been predicted for some models of interfacial roughening transitions in quasi-crystalline 
materials this is not of particular importance to the general fiuctuation theory of wetting 
(and filling) and will not be considered here. The predictions we make for the values of 
critical exponents at filling will assume that in 2D, 1 > C > 1/2. 

2.2. Fluctuation regimes for wetting 

Quite generally, provided C > 0, the critical singularities at wetting are believed to fall 
into three distinct classes depending on the interplay between interfacial wandering and 



the 'direct' infiuence of intermolecular forces ||13| , |Tj, |l^, ^Tj. The latter can be 
modelled by the binding potential 

1^(0 = -^ + ^ + .. .;/>0 (16) 

with a, b effective Hamaker constants and indices q > p > depending on the range of 
the intermolecular forces. The binding potential describes the bare or mean-field wetting 
transition and in order for this to be continuous we require that a = 0, 6 > at the 
(mean- field) phase boundary. Thus a oc (T^lf — T) where T^if is the mean-field wetting 
temperature. The existence of three regimes can be understood semi-quantitatively 
by comparing the bare binding potential with an effective fluctuation contribution 

2^1 



Wfiil) - ^ (17) 



the form of which is suggested by interfacial Hamiltonian models. Beyond mean-fleld 
level we anticipate large scale fluctuation effects with 1^, ~ and thus using (|^) one 
can estimate 

^ ^ l'-'/^ (18) 

^11 
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implying that Wfiil) ~ /^"^^ The competition between Wjiil) and the bare potential 
leads to the following three fluctuation regimes 

• Mean field (MF) regime - If g < 2(1/C — 1) fluctuation effects are negligible, 
K ^ the mean-interface position remains close to the minimum of the 
binding potential. The phase boundary remains a = with f3s = l/{q ~ p). 

• Weak fluctuation (WFL) regime - If g > 2(1/C - 1) but j9 < 2(1/C - 1) the 
repulsion from the wall has an entropic origin but the attraction is still due to 
long-ranged forces. In this regime the phase boundary remains a = with 

ft - WK^. <^'' 

with large-scale interfacial fluctuations It^ ~ 

• Strong fluctuation (SFL) regime - If p > 2(1/C — 1), fluctuations dominate 
leading to a renormalisation of the phase boundary and universal critical behaviour. 
Because the wetting phase boundary no longer occurs at a = one cannot use the 
above heuristic argument to determine the values of the critical exponents. However 
a remarkable feature of 2D wetting is that the values of the SFL regime critical 
exponents can be explicitly related to the wandering exponent ( using very general 
random walk arguments |TB|. For wetting in systems with ( > 1/2 the full set of 



values for the critical exponents is 

«. = 0,/3. = Y3^,z/|| = (20) 

and 

A = ^,V' = ^^ (21) 

1-C 2-C ^ ' 

which are in precise agreement with the explicit results of Ising model [2^ and 



effective interfacial Hamiltonian studies |^ |2^, The values of these critical 

exponents will play a central role in our discussion of 2D filling. 

Finally we point out that in both the WFL and SFL regime, where ~ ^± the 
divergence of the interfacial height, at /i = 0, written in terms of the contact angle, 

W) ~ r^'^ (22) 

is characterised by a universal critical exponent 

4 - ^ = (23) 
2- as 1-C 

which follows directly from the critical exponent relation (|1^) without using the explicit 
values of the critical exponents (p!9D,(pOD. We shall return to this later when we use the 
covariance relations to determine the values of the critical exponents at 2D filling. 
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2.3. Scaling of the PDF and the short- distance expansion 

The position dependence of the equihbrium density profile, p{z), local response functions 
such as the compressibility/susceptibility xi.^) dp{z)/dji and also higher-point 
functions all show scaling behaviour in the WFL and SFL scaling regimes. The scaling 
emerges in the appropriate limits 2;— >oo,t'— i>0,/i^O with zt'^% ht'~^ arbitrary and 
for the profile we anticipate |3( 



p{z) = pi- {pi - p,)E^{zt'('%ht'-^) (24) 

where E.„{x,y) is the scaling function satisfying 3^,(0, y) = 0, 2,^(00, y) = 1 and 
which is distinct in the SFL and WFL regimes. Clearly the scaling of p{z) does not 
include oscillatory structure close to the wall or effects associated with bulk and surface 
criticality but rather reflects the large scale fluctuations of the unbinding liquid- vapour 
interface. For large z but zt'^'' the profile has a characteristic algebraic short- 
distance expansion (SDK) |3^ 

piz) - pi ^ ip, - pi)izt'-^^r (25) 

where, for simplicity we have set h = 0. The SDE critical exponent 7 (referred to as 
6 in earlier work) also describes the behaviour of the local compressibility and pair- 
correlation functions close to the wall. The critical exponent 7 is not independent, 
and can be related to others using standard surface Maxwell conditions and sum-rules. 
Importantly it takes different universal values in the SFL and WFL regimes and can be 
identified as [^, |3( 



^SFL ^ 2(i/C - 1) - 1/(3, (26) 

and 

7^^^ = 2/C-l (27) 

which are valid for arbitrary dimensions. The scaling of the profile is directly related to 
the scaling of the interfacial height PDF, P-jtil), since interfacial fluctuations dominate 
the distribution of matter and we may write 

p{z) = pi- {pi - p,) [ P^{l)dl (28) 
Jo 

where we have assumed the interface separates regions of bulk vapour and liquid density 
and we have omitted the field dependence of the PDF. Throughout this paper we shall 
omit the field dependence whenever the equation containing it is exact within effective 
Hamiltonian theory and not just restricted to the asymptotic critical regime. Since the 
contact angle 9 is an equivalent (possibly non-linear) measure of the temperature-like 
scahng field t' it is possible to write the scaling dependence as P.„{l;9,h) rather than 
PTr{l',t' , h). Moreover for the next few sections we concentrate on behaviour occurring 
at h = and define P^(/;6',0) = Pt,{1;9). Scaling then implies that in the WFL and 
SFL regimes 

P^{1;9) = a9^'A^{al9^') (29) 
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where A^(x) is the scahng function and a is a suitable metric factor having dimensions 
of inverse length. This may be chosen so that the argument of the scaling function 
is equivalent to the length-scales ratio 1/It^{9). Notice that the power- law dependence 
of the contact angle follows from (^3]) and is the same in the WFL and SFL regimes. 
The different fluctuation effects occurring in these regimes are distinguished by the 
appropriate scaling functions A.^^^{x) and A}^^^{x). Both functions are similar at 
large distances where, out of the range of the binding potential, they decay exponentially 
quickly but have quite distinct short- distance expansions Iy-n{x) ~ x'^^^ describing the 
limit ///^ 0. Using the appropriate values for the critical exponents in 2D it can be 
seen that, for both pure and impure systems, the interface makes many more excursions 
to the wall in the SFL regime than in the WFL regime. Hereafter we will only need to 
deal with the properties of the SFL regime. 

The explicit results of effective interfacial Hamiltonian studies are completely 
consistent with the scaling predictions and SDE. For pure systems with just thermal 
disorder the PDF in the SFL regime is particularly simple P, |27[| 

P^{l-e) = 2T.ee-^^'' (30) 

where S denotes the stiffness coefficient of the unbinding interface and may be identified 
with aiv for continuum, fluid-like systems for which the interface is isotropic. For random 
bond disorder the expression for P7r(/; 9) is considerably more complicated but can still 
be calculated analytically using replica trick methods E 



^^^^'^^ ~ ^V^Tn Jo + ^^'^ 

where k is the inverse length-scale associated with the disorder (see later). Note that in 
both these expressions the respective combinations W and W^, together with the SDK's, 
are in agreement with the above scaling theory. 

3. The point tension for 2D wetting in pure and impure systems 

3.1. The Indekeu-Robledo conjecture 

The one ingredient missing in our review of fluctuation effects at 2D wetting is the nature 
of the point tension r measuring the excess free-energy associated with the point of three- 
phase contact between wall- vapour and wall-liquid interfaces [|^, 0, |3^. The 
reason for this, as first pointed out by Abraham, Latremoliere and Upton (ALU) , is 
that beyond mean-field level, fiuctuation effects make the definition of r a rather subtle 
issue. The purpose of this long section is to identify a method of defining r within 
continuum effective interfacial Hamiltonian theory that we can apply to the case of 
wetting with random-bond disorder. To begin recall that within mean-field theory, it is 
straight-forward to to define the point/line tension r by simply subtracting the necessary 
bulk and interfacial contributions from the grand-potential of the heterogeneous wall- 



fiuid interface [19, 20, 31]. Such studies reveal that as T T^ct at h = 0, r contains a 
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singular contribution which we write 

^sing ~ i'^""' (32) 

with a point/hne tension specific heat exponent ai which depends sensitively on the 
range of the forces. A crucial insight into the singularities of the point/line tension was 
made by Indekeu and Robledo |T^ who pointed out that within mean-field theory the 
singularities of r were consistent with the hyper-scaling equation 

Oii = Ois + i^w (33) 

and conjectured that this is generally valid even in the presence of fiuctuation effects. 
The Indekeu-Robledo conjecture is important because it relates the excess free-energy 
of a heterogeneous wall-fluid interface to the properties of a homogeneous wall-vapour 
interface. As we shall show there are very good reasons for regarding this as the precursor 
of a covariance relationship between flUing and critical wetting. Indeed, we shall be able 
to derive a precise relation for the point tension, valid for pure and impure systems, 
which is in perfect agreement with (p3D. Assuming the validity of the hyperscaling 
relation for the point tension in 2D we can identify, in the SFL regime 

2-«, = \^ (34) 

which completes our list of critical wetting exponents. 

Following our earlier discussion of the interfacial height and PDF it is convenient 
for later purposes to measure the singular contribution to the point tension as a function 
of the contact angle rather than t'. We will only consider the point tension for systems 
with short-ranged forces (belonging to the SFL regime) and thus expect 

rsing(^^) ~ e'-'^' (35) 

with 2 — ai given by ([3^ ) In fact, as we shall show, for interfacial models with strictly 
short-ranged (contact) forces there is no ambiguity deflning the point tension itself t{9) 
to be a function of 6, although, of course such models are only well-defined when 6 is 
small. Thus, for pure systems with ( = 1/2, the Indekeu-Robledo hyper-scaling relation 
predicts 2 — = which may either mean that r remains finite at or diverges (or 
vanishes) more slowly than any power law. For random bonds however, the Indekeu- 
Robledo prediction is unambiguous: a; = 3 so that Xsing ~ 9^^. These predictions are 
discussed in detail below. 



3.2. The ALU point tension for pure systems 

ALU point out that, at least for two-dimensional systems, considerable care has to 
be taken in defining the point tension due to the infiuence of large scale interfacial 
wandering which smooths out the point of contact. In mean- field theories, which ignore 
fiuctuation effects, there is no pathology involved in constructing boundary conditions 
which induce a line of contact between a wall- vapour interface and the edge of an infinite 
drop. However, this latter concept becomes ill-defined as soon as fiuctuation effects are 
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introduced since the surface of an infinite drop of liquid has unbounded fluctuations 
in the interfacial height. To overcome this problem ALU propose a fluctuation theory- 
based definition of r, involving a convolution of partition functions. This requires as 
input some appropriate choice for the partition function representing a finite-size liquid 
drop. In this way the thermodynamic limit can be taken yielding a well-defined point 
tension, although the expression obtained depends crucially on the choice of restricted 
partition function used to model the edge of the liquid drop. The most satisfactory 
definition shows the singular behaviour 

rsi„gW~-ln^^ (36) 

close to the wetting transition which is indeed consistent with the Indekeu-Robledo 
exponent relation. The other fluctuation definitions considered by ALU yield point 
tensions that are either non-singular or have a different numerical pre-factor of the 
logarithm (other than unity). However as pointed out by ALU the above divergence is 
appealing since it coincides precisely with the singularity predicted by a heuristic energy- 
entropy balance argument. In an infinite Ising strip of width L lattice spacings and 
with opposing surface fields (which we refer to as a H — strip), it is well understood that 
pseudo-phase coexistence only occurs below the wetting temperature and for sufficiently 



large strip widths ^ ^ 0. This behaviour is characterised by an exponentially 
large correlation length (see their Fig. 3) reflecting the asymptotic degeneracy of the 
lowest two transfer-matrix eigenvalues. Physically this mean that the interface sticks to 
each wall over exponentially large distances which can be estimated by 

^0 

valid for large L, small 9 and 9L oo. Here is an appropriate (non-singular) 
length-scale for measuring distances along the strip which we anticipate is of the order 
of the bulk correlation length. This has to be introduced for dimensional reasons and 
plays no role in determining the divergence of r near wetting. The argument of the 
exponential reflects the free-energy cost of an interface jump from one wall to the other 
with contributions arising from surface free-energies, leading to the T^OL term and two 
point tensions. Exact evaluation of the correlation length in the H — Ising model and 
also in solid-on-solid approximation (valid at low temperatures away from the bulk 
critical point) yields a point tension in precise agreement with result (^) from the ALU 
convolution definition. 



3. 3. Continuum interfacial models of the point tension in pure systems 

The purpose of this subsection is to show that one may also obtain the logarithmic 
singularity of the point tension for pure systems from the asymptotic scaling of the 
PDF as evaluated using a continuum effective interfacial Hamiltonian. The advantage 
of this approach is that it can be readily generalised to systems with random-bond 
disorder, which we will consider next. To begin, we introduce the interfacial model and 
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explicitly evaluate the point tension using the ALU finite-size H — strip identification 
discussed above. 

The general fiuctuation theory of wetting 2D in pure systems, without quenched 
impurities, is based on the interfacial model (see and references therein) 

H[l] = I dxl^^(^^y + hl + W{l)^ (38) 

where l{x) is the local height of the unbinding (liquid- vapour) at position x along the 
wall and W{1), S are the binding potential and stiffness coefficient introduced earlier. 
We emphasise again that we will focus on isotropic bulk fiuid systems and identify 
E with the surface tension ai^. The partition function Zt^{Ii,12; X) of a interface of 
length X with fixed end positions /(O) = l,l{X) = I' is expressed in spectral form using 
continuum transfer-matrix methods [p^ , 

Z.(/,r;X) = 5^C(/')V^„(/)e-^"^ (39) 

n 

where the eigenvalues and eigenfunctions, labelled n = 0, 1,2... satisfy the Schrodinger 
equation 

- + (hi + W{l))Ml) = EnMl) (40) 

Thus in the thermodynamic limit X —>■ oo the singular part to the free-energy is simply 

/sing = Eo (41) 

which, for h = and T < allows us to identify the contact angle 

0-(--^Y" (42) 



Here we have used Young's equation in the small contact limit for which the interfacial 
model is valid. Similarly the normalised interfacial height PDF and parallel correlation 
length follows as 

PAi) = \Mi)\' (43) 

and 

respectively. For future reference we also define the matrix elements 

{m\f{l)\n) = J dirmil)fWn{l) (45) 
which will appear in the transfer matrix theory for the wedge geometry. As discussed by 



Burkhardt [^], the scaling form of the PDF at (|30|) at /i = characteristic of the SFL 
regime, together with the values of the critical exponents quoted earlier (with C, = 1/2), 
readily emerges from the transfer matrix formulism if, instead of the binding potential 
contribution to (^), one imposes the boundary condition on the wave-functions 

V^;(o) = -xMo) (46) 
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where, for pure systems, A oc t' . Specifically, the ground-state energy Eq = — A^/(2E), 
contact angle 6 = A/(S) and mean interfacial height 1^^ = 1/(2S0). 

Next consider the interfacial model of the finite width H — Ising strip at bulk 
coexistence. Since the system has short-ranged forces we can mimic the influence of 
the surface fields through the boundary conditions 

ij'^iO) = -\MO),i^niL) = XML) (47) 

For T„et > T and large L the first two eigenfunctions are 

Ml) oc cosh(v^2S|Eo|)(/ - L/2)) (48) 

and 

?/'i(/) oc sinh(v/(2S|^i|)(/ -L/2)) (49) 
Using the boundary conditions (147|) we find for XL —>■ oo 

Eo = -^{l + ^e~'M...) (50) 

and 

E, = -^{l-4e-'M...) (51) 

Writing these in terms of the contact angle of the semi-infinite geometry we arrive at 
the desired expression for the parallel correlation length 

e,, ~ ^ (52) 

From this we can now extract the desired result for the point tension in the interfacial 
model using the ALU identification (p7|): 

T = -\ne + A (53) 



where A = — \n 2v^S^o may be regarded as an unimportant non-singular contribution, 
the value of which depends on the choice of reference length-scale ^o- In the interfacial 
model (^) with short-ranged forces the only possible choice of length-scale intrinsic 
to the interface is the inverse surface stiffness so that .^o = 1/^ which is directly 
proportional to the bulk correlation length (and recall we have set = 1). With 
this choice of reference length-scale our expression for the point tension in pure systems 
is simply 

r(^) = -ln^-ln2 (54) 
In terms of the parameters S and A this is equivalent to 

r = -ln^-ln2 (55) 

which will be useful when we consider random-bond systems. 

Having derived this result using the ALU identification observe that the logarithmic 
divergence of the point tension for pure systems is also consistent with the behaviour 
of the PDF Pt^II] 9) for wetting at a single wall. To see this recall that, similar to the 
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correlation length (for the H — strip geometry), the scaling form of the PDF Pn-(/; 9) is 
also determined by the surface free-energy and point tension. Ignoring the normalisation 
constraint for the moment, notice that P-n{l\ 0) may be identified as 

P^(/;^) oce-^'« (56) 

where /^(/) denotes the excess free-energy cost of an interfacial configuration 
constrained to be at height I at some arbitrary position along the wall (which we can 
take to be the origin). For asymptotically large distances / ^ Z,r typical interfacial 
configurations determining P.,^{1]6) will have a triangular shape with incident angle 6. 
The free-energy f^{l) can therefore be estimated as 

/^(/) = 2m + 2r (57) 

showing contributions from two point tensions and the surface free-energies. Notice 
this latter term is precisely twice the value of the analogous contribution to the ALU 
correlation length (|37D so we immediately recover the expression for the PDF (|30|) . To 
extract the point tension from Pn{l', 0) we have to bear in mind that unlike the correlation 
length identification ( pT)) the PDF satisfies the addition constraint of normalisation. 
Moreover in 2D one is not capable of distinguishing the normalisation constant from the 
point tension term e"^"^ since the latter term is simply another constant (independent 
of /). Turning this around we observe that the normalisation constant A^(^) appearing 
in the asymptotic scaling form 

P^(/;^) = Ar(0)e-^'(') (58) 

must be related to the (exponential) of the point tension. The ALU identification of r 
through the correlation length is in fact in precise accord with the behaviour of the 
PDF provided we identify 

P^{l-d) = e-^^''-^/i, (59) 

or more simply 

r = -ln(Areo) (60) 

where again is an appropriate intrinsic length-scale introduced for dimensional reasons 
and which plays no role in determining the divergence of r. Using the choice = 1/^ 



appropriate to the interfacial model we recover the ALU identification (|5J). 



3.4- The point tension for critical wetting with random bonds 

We now turn to the evaluation of the point tension for 2D critical wetting with random 
bond disorder. Together with the interfacial model result for pure systems this will be 
crucial in our discussion of the scaling connections between filling and the SFL regime 
of critical wetting. With random bond disorder the interfacial model for 2D wetting is 



written 13, 26 



H[l] = j dx\^ +hl + W{1) + Vr{x, I) \ (61) 
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where the Gaussian random variable V^(x, /) has statistical properties 



Vr{x,l) = Q (62) 



Vr{x, l)Vr{x', I') - Vr{x, I) Vr{x' , I') = A6 {x - x') 5 {I - I') (63) 

where the overbar denotes an average over the quenched disorder with strength A. It 
is convenient to introduce the length-scale n = AS/2 as a measure of the bulk disorder 
which vanishes for the pure (thermal) system. As first shown by Kardar |^B|] the model 
can be studied using the replica trick identification 

= lim ^" ~ ^ (64) 

n^O n 

where may be interpreted as the partition function for n non-interacting interfaces in 
an environment with bulk random bonds. Some details of this calculation are repeated 
below together with the results necessary for the calculation of the point tension and 
later the wedge interfacial height distribution function and free-energy. Performing 
the disorder average introduces interactions described by the many-body Hamiltonian 
(ignoring / independent terms) 



5ik-l,) (65) 



H[{k}] = Jdx{j2 ) +hl + W (/„) j - A ^ . 

so that the interacting n-body partition function for interfaces of length X with 
boundary values /j at x = and /• at x = X has the spectral expansion 

oo 

ZAh,K;X) = (66) 

m=0 

where ifjm^ is the m-th state wave-function with eigenvalue for n-interacting 

interfaces satisfying the Schrodinger equation H^'^^ipm^ = Emipm^ . Again, ignoring 
constant terms, the Hamiltonian operator is 



(67) 

For systems with strictly short-ranged forces, characteristic of the SFL regime, it is 
convenient to adopt the natural generalisation of the boundary condition ( ^6l) which 



reads 13 



or any Ij. Note that the length-scale A is characteristic of the pure wall- fluid 
interface and remains finite at the wetting transition in the presence of random-bonds. 
The groundstate solution to the eigenvalue problem is given by the Bethe- ansat^; 
wavef unction []T3|, ^] 

= C„(A,/t)e-^S'^+'^S-<.l'-'^l (69) 
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with normalisation constant 

where r(a;) is the usual gamma function. From the wave-function one can easily obtain 
the groundstate energy, and by considering the limit of Eq"'^ /n as n — identify 

/■ =-^^ (71) 
as the singular contribution to the free-energy. Thus the contact angle is simply 

(72) 

and note that these expressions identically reproduce the results for the pure system 
when K — i> 0. From the above it is clear that the disorder lowers the wetting transition 
temperature which now occurs at A = k. We shall also need the expression for the mean 
interface height 

= h^i^ - 

which involves the derivative of the psi or digamma function defined by 

Again in the limit k ^ this reproduces the appropriate result = 1/(2A) for pure 
systems. For finite /t however the asymptotic divergence of as A ^ k is different to 
the pure system and 

From the results for the singular contribution to the free energy and divergence of the 
interfacial height we have = and /3s = 2 in agreement with the general expectations 
for the SFL regime with C, = 2/3. The PDF PTr{l; 9) describing the fluctuations of the 
interfacial height in the asymptotic scaling regime is given by (|3T|). 

To evaluate the point tension for random-bonds we use the properties of the 
probability distribution function taking care to extract the relevant quantities at finite 
n before continuing to n = 0. The n— point PDF is the square of the groundstate 
wavefunction which may be written as the ordered product 

n 

pW ({/,.}) = ^^(A, «:) J]e-2(^+("+i-2,»z, ^7g^ 

i=i 

with li < I2 < .. < In- By analogy with the interpretation of the PDF for pure 
systems the coefficient of each Ij term appearing in the exponential may be viewed 
as the surface free-energy cost of constraining the height of the j^^ interface whilst 
the normalisation constant contains the required information about the point tension. 
Using the appropriate replica trick identification we generalise the result (|^) for the 
point tension in the pure system to 

r = - \im -(cl{\K){i^^r - I) (77) 
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where, in an obvious notation is a suitable choice of reference length-scale for the 
random-bond system which plays the same, trivial dimensional role as the length-scale 
^0 for systems with purely thermal disorder. We emphasise that the choice of does 
not influence the asymptotic divergence of the point tensions as T — >^ T^et and only 
contributes towards the non-singular, background term analogous to the constant A 
appearing in ([S^). Thus we find 

A 



r 



which again introduces the digamma function. We now choose the value of so that 
upon taking the limit k — we recover the correct background term for the pure system 
(^). As K ^ the argument of the digamma function diverges and we can use the 
asymptotic large x expansion 

V'(x) ~ Inx - — + .... (79) 
2x 

Note that the necessary logarithmic singularity for the point tension in the pure system 
emerges naturally from the properties of the digamma function. The appropriate choice 
of reference length-scale is therefore = (S)~^ = and is identical to that chosen 
in our earlier discussion of purely thermal disorder. We regard this as a rather pleasing 
feature of the present replica trick definition of r using the PDF. In terms of the length- 
scales A, K and S our expression for the point tension with random-bond disorder is 
therefore 

r = -^f--l)+lnf;^) (80) 



Alternatively for fixed S and k we can eliminate A and rewrite this in terms of the 
contact angle 6 recalling that 6 (x X — k (x t' . This is the form that is most convenient 
for the discussing the connection with 2D filling. Our final result is 



^('"-<t)-'"(I^) 



which should be compared with ( p^ for the pure system. Equations ( pOf ) and ( pT]) are 
the main new results of this section and will play an important role in our discussion of 
2D wedge filling with random-bond disorder. 

We are now in a position to test the validity of the Indekeu-Robledo critical 
exponent relation for the line/point tension. The singularities of the point tension 
occurring as 6' — > are contained within the digamma function which diverges as 
ip ~ — 1/a; as a; —>■ 0. Thus we can identify the singular contribution to the point 
tension 

^.ing ~ ^ (82) 

implying ai = 3 which is in precise agreement with the conjectured exponent relation. 
Note also that as with the pure system the point tension diverges to +oo as T — > T„et 
although the quantitative divergence is much stronger. 
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4. Two dimensional filling in pure and impure systems I: Scaling theory 

Our presentation of fluctuation effects occurring at 2D filling transitions parallels our 
earlier treatment of critical wetting. In turn we will consider (A) the definitions 
of critical exponents and the derivation of exponent relations, (B) a discussion of 
fluctuation regimes from heuristic scaling arguments and (C) the scaling and SDE of the 
density profile and PDF. From these preliminary considerations will emerge a possible 
fluctuation-induced connection with the SFL regime of critical wetting which will be 
precised later using covariance relations. 

4.1. Critical exponents and exponent relations 

A 2D wedge is a 'V shaped substrate formed from the junction of two planar (identical) 
walls that meet at the origin (say) with angles +a and —a measured w.r.t to the 
z = line. Thus the height of the wall above the line is described by a wall-function 
Zw{.x) = cot a I a; I, where the x-axis runs across the wedge. The wedge is considered to be 
in contact with a bulk vapour phase at temperature T and chemical potential /i and is 
supposed to preferentially adsorb the liquid phase along the surface of the walls and, in 
particular, the wedge bottom. Thus the equilibrium density proflle p{z, x) is liquid-Uke 
in the flUed region and will show packing effects very close to the wall although these 
will not be our concern here. Very general macroscopic, thermodynamic arguments 
1^, 1^ indicate that at bulk coexistence /i = /isat the wedge is completely flUed by 
liquid provided the planar contact angle satisfies 6 < a. Thus for the most common 
case where the contact angle decreases with increasing temperature, the transition from 
partial to complete filling occurs at a filling transition temperature Tgu satisfying 

^(Tfiu) = a (83) 

This implies that complete filling precedes complete wetting and also that the 
filling temperature Tgu can be lowered simply by increasing the angle of the wedge. 
On approaching the filling phase boundary at (Tfiu, /isat(^fiii)) the mean height of 
the interface 1^, as measured from the wedge bottom, diverges. The divergence 
is discontinuous and continuous for first-order and second-order (critical) filling 
respectively. Whilst in 3D both types of transition are possible, in 2D filling transitions 
will almost always be continuous. At two-phase coexistence the equilibrium height 
profile leq{x) measured from the z = line, with ly^ = leq{0) is essentially fiat in the 
filled region of the wedge owing to the absence of any macroscopic curvature as dictated 
by the Laplace equation. The lateral extent of the filled region is therefore controlled 
by a correlation length C,x ~ 2/^„ cot a which is trivially related to the interfacial 
height. Critical effects at 2D filling may be viewed as arising from breather-mode- 
like fiuctuations in the interface height which roll the points of contact up and down 
the sides of the wedge thus changing the height and volume of the filled region. A 
similar picture holds for 3D conic filling but is modified in a 3D wedge owing to long- 
wavelength fiuctuations along the system. Similar to wetting the filling transition has 
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two relevant scaling fields which we can write t = (Tgu— T)/7fiii and h = (pz— Pt.)(/^sat~A*) 
respectively. However for filling one has the additional possibility of using the wedge 
angle to control the deviation from the phase boundary. Thus for fixed T close to Tgu 
the combination ^ — a oc t is a linear measure of the temperature-like scaling variable. 
At bulk coexistence = /^sat the divergence of the mid-point interfacial-height and 
roughness are characterised by critical exponents 

(84) 

and we anticipate that in a fluctuation-dominated regime Pw = ^± so that l^ ~ 
Along the filling critical isotherm T = Tgu, h ^ 0, the midpoint height 1^, (and also 
diverges and we introduce the critical exponent 

L ~ h"^"^ (85) 

to characterise this. Two other critical exponents are defined from the singularities of 
the wedge free-energy fw{6, a, h) which, for later purposes, wc have written as a function 
of the variables which highlight the covariance with wetting. Also for /i = wc define 
fyj{6,a) = fu]{6,a,0). At a thermodynamic level the wedge free-energy is defined by 
subtracting from the total grand potential fl the bulk free-energy and the contribution 
from two (infinite) planar walls: 

U{0,a,h) = Q + pV-al^^A (86) 

where awv is the wall-vapour tension for the planar {a = 0) system and A is the total 
surface area exposed to fiuid. By construction the wedge free-energy vanishes in the 
planar limit a = 0. On the other hand we expect that /«,(«, a) is unbounded due to the 
adsorption of a macroscopic amount of liquid. Near the filling transition wc anticipate 
that fw contains a singular contribution that the shows scaling behaviour depending on 
the variables h and t (x {9 — a) only. We write 

which introduces the wedge specific heat exponent Oyj, gap exponent and free- 
energy scaling function Ww{x). Partial derivatives of the wedge free-energy are related 
to thermodynamic observables, similar to the Gibbs adsorption equation for planar 
systems. Firstly, in the free-energy the bulk ordering-field h is conjugate to the total 
2D volume of adsorbed fluid so that 

a /-sing 

^ocC (88) 

Secondly variation of the wedge angle a linearly changes the height and lateral extent 
of the filled region, implying 



oc l^ (89) 
In this way we obtain the exponent relations 

= 2 - + 2/3^ (90) 
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and 

I- a^ = -(3^ (91) 
showing there is only one free critical exponent for 2D filling. 

4-2. Fluctuation regimes for 2D filling 

The classification of fluctuation regimes and also the values of the critical exponents for 
2D filling follow from a rather simple heuristic scaling theory somewhat analogous to the 
Lipowsky-Fisher treatment of critical wetting considered earlier. To begin we consider 
mean-field theory which ignores the fiuctuation effects arising from thermal excitations 
or quenched impurities. As first shown by Rejmer et al interfacial models give 
a very elegant description of filling phenomena at mean-field level. For open wedges 
corresponding to small a (for which tana ~ a) the equilibrium mean-field profile leq{x) 
may be found from minimisation of the effective interfacial free-energy 

F^[l] = j dx\^(^^ +h{l-a\x\) + W{l-a\x\)^ (92) 

which can be justified from analysis of a more general drumhead-like model valid 
for larger a. As mentioned earlier, l{x) denotes the interfacial height relative to the 
2; = line, whilst W{1) is the binding potential appropriate to the planar system. We 
emphasise that the small a approximation is not expected to introduce any peculiarities 
and the critical behaviour predicted by the model (at mean-field level and beyond) 
is believed to be valid for arbitrary wedge angles. The free-energy functional is 
minimised subject to the appropriate boundary conditions that the equilibrium profile 
leq{x) a\x\ + It, as |x| 00. The resulting Euler-Lagrange equation can be integrated 
once to give an explicit equation for the mid-point height (restricting our attention to 
h = 0) i, I, 0, 

- WiL) (93) 

which can solved trivially. Thus for binding potentials of the form (|T6]) the dependence 
on q is unimportant and the mean-field critical behaviour is determined solely by the 
leading order index p. As 6 —>■ a we may expand the above equation 

9 = a + ^lJ + ... (94) 

recalling that a remains finite at the filling transitions since Tfiu < T^et- From this it 
follows immediately that the mean-field value of the height critical exponent at filling is 
= 1/p 0. The absence of any q dependence is a first indication that critical filling 
may be less sensitive to the nature of the intermolecular forces compared to critical 
wetting. 

It is possible to extend this simple mean-field approach to include thermal and 
disorder-induced fiuctuation effects in a heuristic way. The last equation tells us how the 
difference or shift between the contact angle and wedge tilt angle depends on the direct 
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influence of tlie intermolecular forces. In the presence of fluctuation effects arising from 
eitlier tfiermal or random-bond disorder it is natural to suppose that (p^ generalises to 



^ = a + Aap(L) + Aa/i(/^) (95) 

where Aap(/) ~ is the direct angle shift and Aaj/(/) is the shift arising due to 
fluctuation effects which we anticipate takes the form of a ratio of length-scales. Now 
by construction Aafi{l) is not the ratio Iw/^x since this is, essentially, the first term 
in (|95D and is purely geometrical. Instead we write Aa//(Z) oc ^±/C,fi where ^fi is 
an appropriate fluctuation-related length-scale which must be much larger than the 
mid-point height. If we also make the reasonable assumption that this length-scale is 
controlled by the wandering exponent ( then the simplest possible choice given these 
constraints is ^ji ~ lH'^ similar to the relation ~ /^/^ appropriate to a fluctuation- 
dominated (WFL or SFL regime) wetting transition. Consequently if fluctuations 
dominate we anticipate 

Aafiil) ^ l^^^/^ (96) 

similar to the length-scale ratio (HID appearing in the Lipowsky-Fisher analysis. 

This simple, heuristic modification of the mean-field analysis is particularly 
powerful because for filling the phase boundary always remains 6 = a and is not modified 
by fluctuations. It follows that the critical behaviour should fall into two possible classes: 

• Filling mean field (FMF) regime -ifp< 1/(^—1 fluctuation effects are 
negligible, ly^ ^ C,± and the critical exponent (3^, = 1/p is unchanged from its 
mean-field value. 

• Filling fluctuation (FFL) regime -ifp > — 1 there are large-scale 
fluctuations. If ~ C,±, and the critical exponents are universal and determined 
by the wandering exponent. For the divergence of the filling height we predict 

= (97) 

with the values of the other critical exponents following from the relations ( PP] ) and 

m 



At this point, a number of remarks are in order: 

(I) These predictions are in perfect agreement with exact results known from transfer- 
matrix and replica trick studies of interfacial models which find (3^ = 1 and = 2 
for pure {( = 1/2) and impure {( = 2/3) systems with short-ranged forces |^. 
They are also consistent with studies of filling (corner wetting) in square lattice 
Ising models |ll], ^ Moreover for pure systems it is possible to completely 



classify the critical behaviour using the interfacial model and show that the crit- 
icality falls into the above two regimes with a marginal value p = 1 corresponding 
to the FFL /MF borderline. The critical exponent remains = 1 for this marginal 
case. 
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(II) The existence of two fluctuation regimes for filling clearly contrasts with the phe- 
nomenology of critical wetting for which there are three. Also note that the border- 
line between the FFL and FMF regime occurs when p = — l which is different to 
the SFL/WFL and WFL/MF borderlines for critical wetting which happen when 
p = 2(1/C — 1) and q = 2(1/ ( — 1) respectively. The regime in which there is 
universal critical behaviour is broader for filling than for wetting. 

(Ill) The value of the critical exponent = 1/p in the FMF regime is different to the 
value (3s = ^/{q — p) in the MF regime of critical wetting. Therefore when the 
intermolecular forces are sufficiently long-ranged to induce induce mean-field- like 
criticality, there is no apparent connection between filling and wetting. However 
for sufficiently short-ranged forces the predicted value for the critical exponent /J^, 
(pTl), belonging to the FFL regime, is the same as the random- walk result for the 
critical exponent (3s for the critical wetting SFL regime ([20|). This is a first hint that 
there may be some fluctuation- induced connection between the two transitions. 

4-3. Scaling of the PDF and short- distance expansion 

In the FFL regime we anticipate that, within the filled region of the wedge, the density 
profile p{z, x) exhibits universal scaling behaviour related to the scaling of the interfacial 
height PDF. We will focus on the behaviour of the density profile and distribution 
function occurring at the centre of the wedge (x = 0) and define Pw{z) = p{z, 0). Thus, 
for z oo, t 0, h ^ with zt'^™ and ht^^^ arbitrary we expect 



where the scaling function satisfies S^(oo,|/) = 1 and S^(0,?/) = for any y. Notice 
that unlike the case of critical wetting, where one has to distinguish between scaling 
behaviour in the SFL and WFL regimes, the scaling function for fluctuation-dominated 
filling is unique. Associated with the scaling of the profile is a SDK describing the 
algebraic behaviour close to the wall compared to the filling height. At /i = we write, 
analogous to ( ^3]) 



which introduces our final critical exponent 7^, for filling and which is only defined for 
the FFL regime. Similar to SDK exponents for SFL and WFL regime wetting, the value 
of the critical exponent 7^; is not independent and can be related to the other exponents 
defined for filling. To see this, consider that the value of the wall-fluid intermolecular 
potential contains an additional short-ranged contribution of strength Hq localised to 
the bottom of the wedge. In a magnetic (Ising) language this is would correspond to an 
incremental point field at the wedge apex and serves only to introduce a new non-singular 
length-scale proportional to the value of the field. This is useful because differentiation 
of the wedge free-energy w.r.t. Kq yields the value of the density at or near the wedge 
bottom. Now the field Hq is irrelevant, in the renormalisation group sense, and can be 



Pw{z) =pi- {pi - p^)E^{zt^-,ht ^-) 



(98) 



Pw{z)-Pl ^ {py - pi){zt'^^) 



(99) 
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included in the scaling hypothesis for the free-energy by allowing for an additional 
scaling variable h^t^^ which is simply the ratio of relevant length-scales. Differentiation 
of the singular contribution to the wedge free-energy therefore implies that the singular 
contribution to the density at the wedge bottom is simply p^(0) ~ t where we have used 
the exponent relation (pT]). However from the SDE we can also identify p«,(0) ~ t^-^^'^ 
implying that 

7^ = (100) 

which will later prove to be an extremely useful exponent relation. 

The scaling of the profile at filling follows from the scaling of the PDF for the 
mid-point interfacial height, written Pw{l',0,a,h), similar to (p8|). At bulk coexistence, 
h = 0, we simply write 

PUl;0,a,0) = PUl;0,a) (101) 

In the FFL regime we expect that Py^{l;6,a) is characterised by a universal scaling 
function Ayj{x) such that 

P^{1; e, a) = a{e - af^Kifll^O - a)^™) (102) 

where A^(x) is a universal function and the inverse length-scale a is chosen, as 
with A^(a;), so that the argument is simply l/lw Clearly the PDF has the SDE 
A«,(x) ~ x"''"^^. The relationship between the universal scaling function ^^{x) for 
filling and the corresponding function h.^^^{x) for 2D wetting will be central to our 
study. 



5. Interfacial models of 2D filling (II): Exact results and covariance 

5.1. Transfer matrix results 

We begin with the transfer matrix theory of filling in pure systems based on the 
interfacial Hamiltonian 

Hn,[l] = j dxl^i^^ +h{l-a\x\) + W{l-a\x\)^ (103) 

valid for open wedges. It is easiest to assume that the horizontal range is [— X/2,X/2] 
with periodic boundary conditions at the end-points. Note that the model trivially 
recovers the interfacial Hamiltonian for planar wetting H\f\ ( PBD when a = 0. Again we 
emphasise that the assumption of small a is not believed to be in any way important 
as regards the critical behaviour occurring near the filling transition and predictions 
based on the above interfacial model are supported by Ising model studies of filling 
at right-angle corners for different lattice types |jll|. To obtain the partition function 
corresponding to the fluctuation sum over Boltzmann weights it is convenient to make 
the change of variable I = 1 — a\x\ in which case we can re- write the Hamiltonian as 

= 2Sa (le - /(O)) + H[l] (104) 
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where 1^ = l{X/2) denotes the end-point interfacial height (relative to the wall) and 
/(O) = /(O) is the mid-pont height above the bottom of the wedge. Thus the angle a 
enters the partition function only through a local exponential boost factor associated 
with the mid-point height (and end-points). The ensemble average (/(O)) defines the 
equilibrium mid-point height 1^ and from ( |104| ) it is immediately apparent that 



which is a precise version of (|89|). The same relation is also valid in the presence of 



random-bond disorder and will prove useful later. The model can be analysed using 
continuum transfer-matrix methods which yield very general expressions for the wedge 
free-energy and interfacial height PDF, valid for general choices of binding potential. In 
the thermodynamic limit X —>■ oo and in terms of the inner product defined in (^) the 
wedge free-energy follows as 

/^(«,^,/i) = -ln(0|e2^°'|0) (106) 

The probability of finding the interface at height / from the wall at position x along it 
is given by 

p. J . (n|e^^"'|0)C(OV'o(Oe(^°-^"^l-l ...7. 

requiring knowledge of the full transfer matrix spectrum of the planar system. At the 
mid-point {x = 0) however, for which I = Z(0), sjnnmetry considerations simplify the 
expression considerably and 

which only depends on the groundstate properties of the planar problem. This is 
indicative that the mid-point PDF will play a special role in the theory of wedge filling. 
Using these relations it is easy to establish that the filling transition is located at ^ = a 
(and h = 0) in precise accord with the thermodynamic prediction. Moreover the critical 
behaviour falls into two categories in agreement with the heuristic treatment of the 
previous section. For binding potentials with p > 1 the asymptotic criticality is mean- 
field-like with jSyu = 1/p and z/j_ = (1 +p)/2p so that ly^ ^ In the FFL regime 
corresponding to p > 1 the behaviour in the asymptotic scaling regime is universal and 
the same as that found for systems with purely short-ranged forces using the boundary- 
conditions (^). At h = the scaling expressions pertinent to this critical regime are 

p^(l- e-a) = 2S(0 - «)e-2^(^-'^)' (110) 
= In(^-a) -In^ (111) 
corresponding to critical exponents = 1 and 2 — = 0(ln). 
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In the presence of random bond- disorder the generahsation of the interfacial model 
(|6T|) for fining transitions in open wedges is 

j ^\^{^ .W-o^.n.wii-<A.n^v.(.,^ (112) 

and for systems with purely short-ranged forces (and at coexistence /i = 0) the model 
can be solved exactly by extending Kardar's replica trick theory described earlier 
The extension is possible because, similar to ( |104| ) the replicated Hamiltonian is the 
same as the corresponding wetting model apart from a sum over terms 2Sa(/j(0) — If.) 
which may be absorbed into the Bethe ansatz (@]). We omit the details and only quote 
the final results for the infinite wedge. The filling transition occurs at 

A = K + Sa (113) 

which by virtue of (|72D is equivalent to the condition 9 = a. The mean mid-point height 
is given exactly by 

which recovers the pure result ( |109D in the limit k ^ 0. As ^ — a at finite k, the 
interfacial height diverges as 

~ TTTT ^; (115) 

2(A — La — ny 

equivalent to ~ 1/(6' — a)^ and implying that (3^, = 2. The scaling form of the PDF 
describing the asymptotic divergence of the ly^ is 

PUh 0, ") = H^e-'(^-)^-V- r rf, (116) 



TTv^ Jo s + 2l{e-afT.^/K 

Finally the wedge free-energy (at /i = 0) is given exactly by 

Ue,o)^U^^\-4'^\ (117) 



which exactly recovers the pure result (|1 1 1|) as k — 0. As 9 a the free-energy shows 
the singular behaviour 

ping 1^ C^^gN 

implying = 3. For both pure and impure systems the wedge free-energy diverges to 
— oo as ^ ^ a. 



5.2. Covariance laws for filling and wetting 

The above results for 2D fluctuation- dominated filling in pure and impure systems point 
to a remarkable connection with the scaling behaviour occurring for the SFL regime of 
critical wetting. This goes far beyond the identity of the exponents (3^ and Ps suggested 
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Figure 2. Intcrfacial configurations contributing to tlie one-point intcrfacial height 
probability distribution function in two different geometries. Figure a) shows the 
triangular-hke configuration an interface adopts when it is constrained to pass through 
some arbitrary point at height ^ ^ at bulk two-phase co-existence. Figure b) shows 
an intcrfacial configuration in a two dimensional wedge geometry represented in terms 
of the relative height I (see text). Near a filling transition 6* « a and the contributing 
profiles to the respective PDF's, PttH) and Pw{l), are essentially the same. 



by the heuristic scahng theory. For systems with strictly short-ranged forces and at 
bulk coexistence {h = 0) we have established the following covariance relations 

L{9,a) = l^{9-a) (119) 

P^{l-9,a) = P^{l-9-a) (120) 

U0,a) = r{9)-r{9-a) (121) 

The final relation between the wedge free-energy and the point tension has not been 
reported before and is one of the central new results of our paper. These 'laws' are 
also valid in the asymptotic critical region, 9 a, even in the presence of long-ranged 
forces provided the filling transition belongs to the FFL regime. We emphasise that 
the connection between filling and the critical wetting SFL regime is all the more 
remarkable because the FFL regime is broader. For example, recall that with purely 
thermal disorder the FFL corresponds to binding potentials with p > 1 whilst the critical 
wetting SFL regime corresponds to p > 2. Thus for model systems with 1 < p < 2 the 
filling transition precisely mimics the properties of the SFL regime even though the 
wetting transition for the corresponding planar system belongs to the WFL regime. It 
is in this sense that the wedge geometry effectively turns off the influence of the long- 
ranged forces. We conjecture that the above covariance relations laws connecting filling 
and wetting are generally true in 2D provided the wandering exponent 1 > C ^ 1/2- 

The fiuctuation-induced covariance between filling and wetting has a simple 
geometric interpretation. In figure ^ are shown typical intcrfacial configurations 
contributing to the PDF for two different geometries. On the LHS is shown the typical 
triangular configuration an interface adopts at a planar wall when it is constrained to 
pass through a point at height / far in excess of the mean intcrfacial height On the 
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RHS is the typical configuration for an interface in a 2D wedge geometry plotted in 
terms of the relative height / = l{x) — a\x\. Close to the filling transition 9 a and 
consequently the typical interfacial fluctuations contributing to the Pjil) and P7r(/) in 
the different geometries are essentially the same. Similar remarks also apply in 3D for 
the cone geometry 0. 

The covariance relations are extremely restrictive and contain a great deal of 
information about the allowed values of the critical exponents at 2D filling and wetting. 
Indeed, it is worthwhile developing the consequences of these relations assuming only 
their validity together with the critical exponent relations derived earlier from standard 
scaling theory. Firstly, from the first relation (|119|) , it necessarily follows that 

= 4 = (122) 

in agreement with the heuristic scaling theory. It is important to realise that this 
identification does not depend on the specific values of the critical exponents pertinent 
to the SFL regime since the divergence of lT:{d) as 6' ^ is determined by the critical 
exponent (3s (0) rather than by /3<j. This is because the covariance relations are 
expressed in terms of the angles 9 and a rather than the scaling fields t and t' . Using the 
derived exponent relations for filling we can now deduce the values of the other critical 
exponents in the FFL regime 

1 — 2C 1 

2 - = ^, = J— ^, = C (123) 

which are all universal, determined only by C,. The second covariance law for the PDF's 
contains even more information. In terms of the scaling functions this reads 

A^(x) = Af^lx) = A(x) (124) 

which clearly indicates that the connection with the SFL regime is fundamental and not 
a merely fortuitous coincidence of critical exponent values. The identity of the scaling 
functions now has a truly remarkable consequence. To see this note that it necessarily 
follows that the SDK exponents 7^ and 7'^^'^ have the same value. Recalling the general 
critical exponent relations (|2^ ) , ( P^) for these leads to the identification 

1//?^ = 2(1/C - 1) - l/(3s (125) 

and using the above value for we find 

= (126) 

which re-derives the random-walk predictions for the critical wetting SFL regime. Thus 
the covariance relations severely restrict the allowed values of the critical exponents for 
both 2D filling and 2D critical wetting. 

Next we turn our attention to the free-energy covariance law. Taking the derivative 
of ( |121D w.r.t. a we find 
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and setting a = we arrive at a novel result relating the planar interfacial height to the 
point tension: 

W) = (128) 

valid for interfacial models with short-ranged forces (or in the asymptotic SFL critical 
regime). This relation has a number of consequences. First, equating the power-law 
critical singularities on either side yields 

2(3s _ 2(2 - ai) 
2 — as 2 — as 

which reduces to 



1 (129) 



ai = as + i>\\ (130) 

thus deriving the conjectured Indekeu-Robledo exponent relation. We emphasise that 
in this manipulation we have only used the general exponent relation (|n|) without 
having to introduce any results specific to d = 2. This is strongly suggestive that a 
generalisation of (|128D, with possibly different numerical pre-factors, may well exist for 



SFL regime critical wetting in higher dimensional systems. Secondly, beyond simple 
power-law singularities we can now see that the logarithmic divergence of the point 
tension for purely thermal systems found by ALU is, in fact, necessary in order that the 
interfacial height diverges as l-,^ ~ 9~^. 

The values of all the other critical exponents for filling and wetting now follow from 
standard exponent relations. The connection between them can be summarised by 

(3^ = Ps, aw = au = i^w (131) 

where the LHS and RHS refer to the FFL and SFL regimes respectively. 

Finally, for completeness, we remark that for pure systems it has been shown 
that for the marginal case p = 1, q = 2 corresponding to the FFL/FMF boundary, 
the covariance laws for the interfacial height and PDF relate the behaviour at filling 
to the WFL/MF regime of wetting. The wedge free-energy can be easily calculated for 
this case and is similar to ( |111| ) but has a numerical pre-factor c = 2 -|- (1 + 8S6)^/^ 



in front of each logarithm. However we do not discuss the possible connection with 
the point tension because we are not confident that for such long-ranged forces r is a 
well-defined quantity. Whilst an expression for r can be found for such systems using 
the PDF identification, and is in accord with the covariance law, we have not been 
able to extract r using another, independent, method. We feel such a check is necessary 
since, as shown by ALU, even for short-ranged forces, the convolution definition of r can 
lead to different results. Moreover for p = 1 it is not obvious that r can be extracted 
using a generalisation of the ALU correlation length identification. It may even be 
that the covariance relation between the wedge-free energy and the point tension can 
be forwarded as a suitable definition of the point tension for interfacial models with 
this marginal interaction. This would certainly be consistent with the Indekeu-Robledo 
conjecture for the point tension singularity. 
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6. Scaling and covariance for the local compressibility 

Our treatment so far has concentrated on critical singularities occurring at bulk two- 
phase coexistence /i = 0. Given the precise connection between 2D filling and wetting 
occurring at coexistence it is natural to enquire whether this extends to quantities 
defined for /i > 0. Away from two-phase coexistence, however, any possible relation 
between 2D filling and wetting is certainly subtler than that occurring for h = because 
the pertinent gap exponents and A are different. Thus the divergences of 1^ ~ h'"^^ 
and ~ along the respective filling and wetting critical isotherms are quite different 
and preclude a law of type ( |119|) . Similarly there can be no simple generalisation of the 
free-energy relation ( |121| ) because the point tension is only defined for h = 0. This 



suggests that we first look for covariance relations between response functions for FFL 
filling and SFL critical wetting describing infinitesimal deviations from bulk coexistence. 
Such relations, should they exist, will also be notable because response functions are 
generally related to integrals over two-point functions which would suggest that these 
too satisfy covariance relations. 

Before we calculate the scaling expressions for PDF and local compressibility for 
2D filling we recall some pertinent results known for critical wetting. 

6.1. SFL critical wetting 

Differentiating the scaling ansatz for the profile p{z) w.r.t. h immediately implies that in 
the scaling limit of the SFL regime, and up to an unimportant non-universal pre-factor 
D, the local compressibility Xn = dp{z)/dh evaluated at bulk-coexistence {h = 0) has 
the form 



xA^; e) = De-'^/^'-^'^x^^\zMe)) (132) 

where X.^{x) is a scaling function describing the universal position dependence. This 
is independent of the range of the forces and is specified by the dimension and type of 
disorder only. Again we emphasise this is valid in the asymptotic scaling limit 6^0, 
z —y oo with z/Itt arbitrary. The SDE is controlled by the same exponent as the density 
profile so that 

X^^\x) ~ x^'^' (133) 

as X — s> 0. By adopting the convention that the pre-factor is unity, the scale of X.^ 
is fixed and we can regard the scaling function as universal. In effective Hamiltonian 
theory the behaviour of the compressibility is directly related to that of the PDF since 
from (|28|) we have 

dP (I) 

XM = {PI-P.) I dl^^ (134) 

The scaling form of the PDF, density profile and local compressibility emerges naturally 
from the interfacial model if we use the same boundary condition (^61) but retain the hi 
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term in the Hamiltonian. As shown by several authors, for /i > the ground-state wave 
function is an Airy function [R^, HO, El| implying 



P^- 9, h) oc Ai^ (h^lH - e^h-^/^W^{he-^)^ (135) 

where, for the sake of clarity, we have dropped non-universal metric factors and written 
the field dependence in terms oi 9 ~ t' . Taking into account the h dependence 
coming from the free-energy scaling function Wt^{x) and the normalisation constant, 
it is straightforward to show that x-k{.z; 6) scales according to the prediction (|132|) with 
a universal scaling function 

Xf^^(x) = {x + \x^)e-^ (136) 

Notice that the SDE behaviour of this function is in accord with the general 
requirement (|133|) although a different power-law determines the algebraic correction 
to the asymptotic exponential decay. 

6.2. FFL filling 

In zero field the mid-point local compressibihty, Xw{.z) = dpw{z)/dh , should also show 
scaling behaviour analogous to the behaviour occurring at SFL regime critical wetting. 
Following our treatment above it follows from the profile equation (^Sp that ioi z oo, 
^ oo with fixed z/l^, 

XU^; 9, a) = b{9 - X^{z/U9, a)) (137) 

where Xyj{x) is the appropriate scaling function whose SDE is described by the critical 
exponent 7^. Thus we expect 

X^{x) ~ x^- (138) 

and by again adopting the convention that the critical amplitude is exactly unity we 
can fix the scale of the universal function X^(a;). 

The behaviour of the mid-point PDF at filling can be easily calculated using the 
transfer matrix result (|108|) and after a little algebra we obtain 



P^il- 9-a,h)(x e-2S(e-a)/-MV" (139) 

which is considerably simpler then the planar result. Note that the h dependence of 
this scaling function has a simple geometrical meaning since the term I'^/a corresponds 
precisely to the area of the filled region of the 2D wedge. From ( |139| ) it is easy to check 
our earlier prediction that, along the critical filling isotherm (T = Tgu,/;, — > 0), the 
interfacial height diverges as 1^ ~ h~^. We find 

L ~ i^f (140) 

in perfect agreement with the expected behaviour for purely thermal disorder = 1/2). 
It is also immediate from ( |139| ) that the roughness ~ 1^ along the critical isotherm. 
We mention in passing here that this behaviour contrasts with the singularities occurring 
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at complete filling corresponding to /i for T > Tgu. For this transition it is apparent 
from ( |139| ) that the critical behaviour is not fluctuation-dominated and Z^, ~ whilst 
~ h"^/"^ . The divergence of 1^ agrees with predictions based solely on thermodynamic 
and mean-field arguments which remain valid because the fluctuation effects at complete 
filling are small even for the present two-dimensional system with short-ranged forces. 
Returning to the local compressibility we can readily calculate the desired expression 
for the zero- field local compressibility using the PDF (|139D . Taking care to account for 
the h dependence of the normalisation factor we find that the local compressibility is 
precisely of the form ( |137|) with = 2 and a scaling function 

X^{x) = (x+^x2)e-^ (141) 
which is identical to that derived for SFL regime critical wetting. Thus we write 

X^{x) = X^/\x) = X{x) (142) 

In terms of the full angle dependence the scaling local compressibilities satisfy the 
simple covariance relationship 

Xw{z;0,a) ^ 0-a 

valid in the asymptotic critical regime 6 —>■ a. The identity of the local compressibility 
scaling functions for filling and SFL wetting is the main result of this section. Whilst 
we have only demonstrated this for pure systems we expect this is also valid for systems 
filling and wetting in other 2D systems provided that ( > 1/2. Support for this 
conjecture comes from the derived values of the critical exponent relations. In particular 
in the limit z/lir —>■ ^ and z/ly^ —>■ the SDK's for the local compressibilities are identical 
and hence their ratio must be independent of z. Moreover the values of the exponents 
and 2 A/ {2 — as) are always such that the ratio is proportional to 6* — a independent 
of C- It therefore seems highly likely that the covariance relation ( |143| ) is also valid for 
impure systems and can be checked in future studies. 



7. Wedge Covariance for pure systems revisited 

Our discussion of the implications of the covariance relations in §5 was simply based on 
comparison with the predictions of standard scaling theory for critical exponent relations 
and the short-distance expansions. With a little more input concerning the properties 
of the interfacial Hamiltonian model we can also use the wedge-covariance relations to 
re-derive the results of the transfer matrix studies without explicit calculation. The 
following discussion is restricted to thermal disorder although a generalisation to the 
impure systems may well be possible using the replica trick method. 

In terms of the relative interfacial height /— we can write the wedge Hamiltonian 
as a perturbation from the planar model 

hJ] = H\1]+6HJ] (144) 
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where 

5H^] = -2Sa j dx5{x)l{x) (145) 

and we have ignored /-independent terms and the boundary conditions concerning the 
end-points /g which do not matter in the thermodynamic hmit. Because the planar 
interfacial Hamihonian and the wedge perturbation are both local it immediately follows 
that, up to a trivial normalisation constant, the mid-point height PDF for pure systems 
satisfies 

P^il) oc P^(/)e2^"' (146) 
which is equivalent to ( |108| ). Wedge covariance of the PDF at bulk coexistence implies 

P^{l;e -a) (X P^{l;e)e^^''^ (147) 
and setting a = 6' it follows that 

P^(/;0) oc P^(/;0)e-22^' (148) 

where the position dependence of Pn{l', 0) is determined solely by the SDK. Using the 
known critical exponent relation for ^-^^^ we conclude that within the SFL regime the 
PDF is necessarily of the form 

P^(/;^)oc/2/f'i/^^-=^e-2^^' (149) 

From this it follows that $s = ^ and since $s = C/i^ ~ only value of ( consistent 
with the wedge covariance hypothesis in pure systems is 

C = I (150) 

implying that /5s = 1. Similarly the universal scaling function for the PDF in the SFL 
must be simply 

A(x) = e"^ (151) 

Now consider the PDF off coexistence and note that the exponential boost factor 
does not depend on the bulk field h. It follows that in the asymptotic critical regime, 
the scaling functions for filling and wetting satisfy the simple quotient relation 

P^{l;e,a,h) ^ P^{l;a,h) 

PUl;0,a,0)^ P^il;a,0) ^ ^ 

where the constant of proportionality is trivially determined from the normalisation 
conditions on each PDF. Keeping the value of ( arbitrary for the moment it follows that 
to first order in h the quotient has the expansion 

where is a pure number and {l^^'^) denotes the appropriate moment of the interfacial 
height evaluated at h = 0. The various terms in this expansion arise for the following 
reasons: (1) the quotient must tend to unity as /i ^ 0; (2) the expansion must be linear 
in h because PT,{l;a,h) is a non-singular function of /i for a > 0; (3) the power-law 
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dependence /^/^ is necessary since along the filling critical isotherm [a = 6, h ^ 
scaling demands that ~ /i"^; (4) the dependence on a follows from conditions (2), 
(3) given the anticipated scaling form of PDF PTt{l\ 6, h); (5) the additive term involving 
(/^/'') arises from the normalisation condition on each PDF. 

From ( |153[ ) we immediately observe that the universal scaling functions A(x) and 
X{x) for FFL filling and SFL wetting must be related according to 

X'(x) oc A(x)(l -c^x^/^) (154) 

where 

/„°°.VCA(,)rf, ^^^^^ 

and the constant of proportionality in (|154|) is trivially fixed by the condition that 
X(x) ~ x^^'^~^ for small x. At this point we can substitute the appropriate expressions 
A(x) = and ( = 1/2 for pure systems and integrate to find 

X(x) = (x + ^x2)e-^ (156) 



in agreement with the transfer matrix calculation. The differential equation ( p.54| ) 
is also an effective method of calculating the covariant scaling function for the local 
compressibility at the FFL/FMF (or equivalently WFL/MF) borderline with p = l,q = 
2 for which the PDF function A(x) has a non-trivial SDK Further work is required 
to see if similar approaches can determine the allowed values of (, A(x) and X(x) for 
disordered systems without using the full transfer matrix formulism. 



8. Discussion 



In this paper we have investigated fluctuation effects occurring at 2D flUing and the 
fundamental connection with the SFL regime of critical wetting. Our study has revealed 
that in addition to the known covariance of the interfacial height and PDF f^, 
the wedge free-energy and mid-point local susceptibility are also related to the point 
tension and local susceptibility at wetting. These relations are extremely restrictive and 
determine the allowed values of the critical exponents at FFL flUing and SFL wetting 
without further assumptions other than those of standard scaling theory. Moreover 
if wedge covariance is combined with knowledge of how the wedge Hamiltonian is 
perturbed from the planar interfacial model, then very speciflc predictions for the value 
of the wandering exponent, critical exponents and scaling functions can be obtained. 
Thus wedge covariance appears to play a similar role to the principle of conformal 
invariance for 2D bulk critical phenomena in that it yields predictions over and above 
those of scaling and scale invariance. Wedge covariance also bring new insights into 
the nature of the SFL regime wetting transition itself. In particular we have shown 
that the covariance relation for the wedge free-energy provides a means of deriving the 
conjectured Indekeu-Robledo relation for the critical singularity of the point tension 
and also explains its logarithmic divergence for pure systems. Other aspects of our 
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work that have not been presented before include the derivation of an expression for the 
point tension for impure systems and the development of a scaling theory and derivation 
of critical exponent relations for 2D filling transitions. 

It is hoped that the present work motivates the further study of filling transitions 
in both two and three dimensional, pure and impure systems. Within the framework 
of the present, small a interfacial model, future work could include discussion of two- 
point functions. These may well exhibit some covariance properties as suggested by 
the behaviour of the local compressibility. Non-linear functional renormalisation group 
analyses, both approximate and exact decimation |]13| type, would also shed light on 
the likely expanded space of the SFL regime fixed point. Studies of fluctuation effects 
at filling using other models such as the full drumhead interfacial Hamiltonian or lattice 
Ising model, for which only limited results are known would also be very welcome. 
It would also be interesting to see if the formal statistical mechanical theory of fluids at 



interfaces, with its powerful sum-rule and correlation function hierarchies 42, 43 



can be applied to the problem. Perturbative expansions of the many-body Hamiltonian 
analogous to (|144|) and ( |145| ) may prove revealing. Staying within the framework of 
effective Hamiltonian theory it would be highly informative if the nature of both 3D 
wedge and cone filling were understood beyond the case of purely thermal disorder. 
It is likely that this would also shed more light on the nature of the SFL regime in 
higher dimensions about which very little is currently known. Nevertheless the simple 
heuristic picture of how the filling of a 2D wedge manages to precisely mimic the 
behaviour of the SFL regime, does generalise rather naturally to the cone geometry 
which similarly enforces the same qualitative type of conic interfacial configuration that 
determines large deviations in the one-point function at wetting. It is tempting to 
speculate that the Indekeu-Robledo critical exponent relation may also be intimately 
tied to a possible covariance for three dimensional wedge free-energy although much 
further work is required to quantify this. 

AOP and AJW are very grateful to Dr. Carlos Rascon for discussions and 
collaboration on related work. AJW and MG thank the EPSRC for financial support. 
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